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ABSTRACT 

The purpose of this presentation is to show that it is possible to use 
nonsmooth optimization algorithms to design both closed-loop finite 
dimensional compensators and open-loop optimal controls for flexible 
structures modeled by partial differential equations. 

An important feature of our approach is that it does not require modal 
decomposition and hence is immune to instabilities caused by spillover 
effects. Furthermore, it can be used to design control systems for struc- 
tures that are modeled by mixed systems of coupled ordinary and partial 
differential equations. 


260 



SCOLE 87/3 


DESIGN OF STABILIZING FEEDBACK-SYSTEM COMPENSATORS 

The optimization-based design of finite dimensional compensators for 
systems modeled by coupled systems of ordinary and partial differential 
equations is made possible by a generalization of the following necessary 
and sufficient stability test for linear systems described by ordinary 
differential equations. 


THE DYNAMICAL SYSTEM 

Consider a parametrized, linear, time-invariant, interconnected, finite 
dimensional dynamical system, X(p), described by a set of state equations: 

x(t) = A(p)xj(t) + B(p)u(t), 
y(t) = C(p)x(t) + D(p)u(t), 

We shall denote the characteristic polynomial of £(p) by X( S >P) and 
assume that the coefficients of x(s,p) are continuously differentiable in p. 
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S-STABILITY 

When, it is desired to ensure not only exponential stability of a closed 
loop system, but also to exercise some control over the location of its 
poles, it is convenient to make use of the following definition of S- 
stability. 

Definition (S -stability): Consider a linear, time-invariant, finite dimen- 

sional dynamical system £ of the form (1). Let S be an open unbounded 
subset of C which is symmetrical with respect to the real axis, and such 
that S c 3 C + , where S c is the complement of S and C + is the closed right 
half of the complex plane. 

We say that the system L is S-stable if all the zeros of its charac- 
teristic polynomial are in S. ■ 
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A MODIFIED NYQUIST STABILITY CRITERION 

Theorem : Let Sc C be as specified in the Definition and let Be C be 

any simply connected set satisfying (0,0) £ B. Suppose that 
D(s,q) e C[s] is a parametrized polynomial of degree N, whose 
coefficients depend on the parameter vector q g R D in such a way that 
for every x(s) e P N satisfying Z[x(s)]cS, there exists a q x e R n ° such 
that 

(i) Z[D(s,q x )]cS, (2a) 

(ii) x(s)/D(s,q x ) g ^ s G ( 2b ) 

Then, given a polynomial x( s ) G Pn» Z[x(s)]cS if and only if there exists 
a q x e R n ° such that (2a, b) hold. ■ 


263 



SCOLE 87/6 


PROOF OF MODIFIED NYQUIST STABILITY CRITERION 

(= = =>) Suppose that Z[x(s)]c:S. Then, by assumption, there exists a 
q x e R nu such that (2a), (2b) hold. 

(< = = = ) Next, suppose that (2a), (2b) hold. Then, because B is a sim- 
ply connected set which does not contain the origin, the locus traced out 

in the complex plane by x( s VD(s,q x ), ^ or s € < * oes n0t enc * rc ^ e 
origin. It now follows from (2a) and the Argument Principle that 

Z[ X (s)]cS. ■ 

Comment : It is clear from the Theorem that an acceptable parametri- 

zation of the polynomial D(s,q) depends on the shape of the set S and the 
choice of the set B. A further requirement is imposed by semi-infinite 
optimization: the parametrization must be such that it is easy to ensure 
that the zeros of D(s,q) are in S. ■ 
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OPTIMIZATION-BASED CONTROL SYSTEM DESIGN 



DESIGN VECTOR: x = [ x\ x 2 ,..., x 8 ] . 
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DESIGN CRITERIA 

1. The feedback system must be exponentially stable. 

2. The system should have a good step input response. 

3. There should be little interaction between channels. 

4. Plant should not be saturated by command input effects. 

5. System should have high output disturbance rejection. 
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MODIFIED NYQUIST STABILITY CONSTRAINT 
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CHANNEL INTERACTION CONSTRAINT 



Magnitude of H y : r i(jco,X) 



Magnitude of H ^(jo^x) 


COMMAND INPUT SATURATION CONSTRAINT 



Singular Values of H U2r (jCD,x) 
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OUTPUT DISTURBANCE SUPPRESSION CONSTRAINT 
Must accept some disturbance amplification outside operating 
bandwidth: 

o[ H yd (jco,x) < 1.05, V © e [1,1000] 

COST: OUTPUT DISTURBANCE SUPPRESSION 

Suppress disturbance effects inside operating bandwidth: 

f(x) = max o[H yd (jco,x) 

©e [0.001,1] 



fr*qu«Bcy (lad/Stc) 


Singular Values of H y( j(j®» x ) 
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INTEGRATED STRUCTURE-CONTROL-SYSTEM DESIGN 


BEAM SHAPE AT TIME T 



Vibrating Beam 
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DYNAMICS 

• GENERAL MODEL: Euler-Bernoulli Model, Kelvin-Voigt or 

Proportional Damping, Coupled Axial and Flexural Linear PDE’s. 

• Control Forces F'(t), Actuator Positions a 1 , Sensor Positions s'. 

• SIMPLIFIED MODEL: Decoupled Motion Formulation: 

n i 

mu„(t , x) + CIu lxxxx + EIu xxxx (t , x) = £ b,(x - 3i ) F'(t) . 

i=l 

1 1 

y'(t) = / c ; (C - s i )u(t , C)dC or y'(t) = J d s (C - s')u(t , C)d£ . 

0 0 

BOUNDARY CONDITIONS 

u(t , 0) = 0 , u x (t , 0) = 0, Ju Ux (t , 1) + CIu lxx (t , 1) + EIu xx (t , 1) = 0 , 
Mu„(t , 1) - CIu, xxx (t , 1) - EIu xxx (t , 1) = 0. 
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DESIGN CRITERIA 

1 Th e feedback system must be exponentially stable. 

2. Control system compensator should be finite dimensional. 

3. Actuators should not be saturated by command input effects. 

4. System should have high mechanical disturbance rejection. 

5. Average power use should be low. 

6. Structure weight should be low. 

7. Structure should remain in elastic range. 
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DESIGN VARIABLES 

• CONTROL SYSTEM COMPENSATOR 

(i) Coefficients of compensator differential equation. 

• STRUCTURE 

(i) Positions of actuators and sensors. 

(ii) Parameters of damping devices. 

(iii) Parameters of composite materials. 

(iv) Parameters determining shape of structure. 
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PRELIMINARY RESULTS 

1. The control system can be stabilized using a finite dimensional 
proportional-plus-integral controller which ensures good distur- 
bance rejection. The use of our modified Nyquist stability cri- 
terion in the design of a stabilizing controller requires only 
evaluations of the system frequency response. Since the frequency 
response at a given frequency can be computed in some cases by 
formula and in the more general cases by solving two-point linear 
boundary value problems, there is no need for modal decomposition 
and hence there are no spillover effects. As in the finite dimensional 
case, time and frequency domain constraints can be treated simul- 
taneously and, in an integrated design approach structural parameters 
and constraints can also be introduced into the optimization problem. 

2. If a sequential design approach is used, an infinite dimensional com- 
pensator can be designed using an H°° frequency domain constraint 
formulation which results in a convex optimization problem and 
automatically ensures exponential stability with stability margin. 

3. An infinite dimensional controller designed as above can be 
approximated by a finite dimensional controller without spillover 
effects. 

4. A special semi-infinite optimization algorithm has been developed 
which is highly effective for design with H°° frequency domain 
design constraints. 
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THE DYNAMICAL SYSTEM 

Hollow aluminum tube: one meter long, 2.0 cm diameter, 1.6 mm 
thick. Attached mass weighs 1 kg. We assume that motor torque u(t) can 
be directly controlled. 

Standard Euler-Bemoulli tube equations with Kelvin-Voigt visco- 
elastic damping: 

mw„(t , x) + CI\v txxxx (t , x) + EIvv xxxx (t , x) - mfi 2 (t)w(t , x) 

(la) 

= - mu(t)x, x g [0,1] 

with boundary conditions: 

w(t,0) = 0, w x (t,0) = 0, CIw txx (t, 1) + EIw xx (t, 1) = 0. (lb) 

M(Q 2 (t)w(t , 1) - w tt (t, 1) - u(t)) + CIw txxx (t, 1) + EI\v xxx (t, 1) = Ole) 

where w(t,x) is displacement of tube from shadow tube (which remains 
undeformed during the motion), u(t) is motor torque, and Q(t) rad/sed is 
angular velocity. Above: m = .2815 kg/m, C = 6.89x1 0 7 pascals/sec., 
E = 6.89x1 0 9 pascals, I = 1.005 x 10 -8 m 4 , The tube is very lightly 
damped (0.1 per cent ). 
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THREE OPTIMAL SLEWING PROBLEMS 


Minimize the time required to rotate the tube 45 degrees, from rest to 
rest, subject to the torque not exceeding 5 newton-meters. 


Minimize the total energy required to rotate the tube 45 degrees, from 
rest to rest, subject to the torque not exceeding 5 newton-meters and 
the maneuver time not exceeding a given bound. 


Minimize the time required to rotate the tube 45 degrees, from rest to 
rest, subject to the torque not exceeding 5 newton-meters and an 
upper bound on the potential energy due to deformation of the tube 
throughout the entire maneuver. 
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THE DYNAMICAL SYSTEM 

MATHEMATICAL FORMULATION OF THE THREE PROBLEMS 

• To avoid technical problems associated with variable intervals and 
problems due to discretization, augment dynamics by one state variable 
and introduce scale factor T > 0 so that problem becomes defined on nor- 
malized time interval [0,1], with T also equal to final time. 

• Tube is at rest when the total energy = energy due to rigid body 
motion + energy due to vibration and deformation = 0. 

(i) To ensure a slewing motion of 45°, we define 

gVT) 4 (e-n/4) 2 ( 2 ) 

(ii) Rigid body energy at final time is proportional to the square of the 
angular velocity. 

g 2 (u,T) 4 O(T) 2 . (3) 
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(iii) Kinetic energy due to tube vibration at normalized time 1 is 

1 

g 3 (t,u) 4 —J w t (l ,x) 2 dx. (4) 

2 0 

(iv) Potential energy due to tube deformation at normalized time 1 is 

g 4 (l,u) = w xx (l,x) 2 dx. (5) 

z 0 

• Potential energy due to deformation of the tube at normalized time t: 

P(t,u) 4 Jp-J w xx (t,x) 2 dx. (6) 

1 0 

(v) To limit tube deformation for all t e [0,1] we define 

T 

g 5 (u,T) = J [max{ P(t,u) — f(t) , 0 }] 2 (7. 

0 

(vi) To ensure slewing time does not exceed T f seconds, we define 

g 6 (u , T) 4 T-Tf. (8 
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FINAL PROBLEM FORM 

P]: Te iRr! n u e G {g0(U ’ T)l8i(U ’ T)_e - 0 -j e {1 ,2,3,4 ) }, 
where g°(u,T) 4 T, R + 4 { y e ]R I y > 0 } and 

G = (ue L„[0, 1] I lu(t)l < 5, Ie [0,1] }. 

?2: T e nTue G { S° (u » T ) 1 g J (u, T ) - £ < 0, j e {1,2, 3, 4, 6 ) ), 

1 

where g°(u,T) ^ j ||u(t)|J 2 dt. 

0 


min 

TelR+.ue 


G * g°(u,T) I g>(u,T) -e SO, j e { 1,2, 3, 4, 5 } }, 


where g°(u , T) = T. 


• All g> are continuously differentiable in LJX), 1]. 
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THE DYNAMICAL SYSTEM 
COMPUTATIONAL RESULTS 

IMPLEMENTATION. Because we cannot solve the system PDEs 
exactly, we cannot evaluate g*(u,T) or Vg-*(u,T) exactly. Furthermore, 
since u is an infinite dimensional design vector, it can only be entered 
into a computer in discretized form. We use an implementable algorithm 
which adjusts integration precision and control discretization adaptively. 
To discretize the PDE in space, we use the finite element method. Since 
the PDE is fourth order in space, it is necessary to use elements of at least 
second order. We have chosen Hermite splines as basis elements. The 
input u e G is discretized in time and Newmark’s method is applied to 
evaluate the resulting system of ordinary differential equations. 

LINEARIZATION. The results presented are for the case in which 
the Q 2 (t) terms are neglected in equation (1). Similar results have been 
obtained by performing experiments when the Q 2 (t) terms are included. 
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TIP DISPLACEMENT FOR MINIMUM-TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY 


Meters 
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DEVIATION FROM SHADOW BEAM FOR MINIMUM-TIME 
PROBLEM WITH TORQUE CONSTRAINTS ONLY 


w(l,x) (meters) 
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CONSTRAINT VIOLATION FOR MINIMUM-TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY: 
DISCRETIZATION EFFECTS 


PSI 
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OPTIMAL TORQUE 

FOR MINIMUM-CONTROL-ENERGY PROBLEM 
WITH TORQUE CONSTRAINTS AND FINAL TIME < 0.8 SEC. 


Control 
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OPTIMAL TORQUE 

FOR MINIMUM-CONTROL-ENERGY PROBLEM 
WITH TORQUE CONSTRAINTS AND FINAL TIME < 1.0 SEC. 
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POTENTIAL ENERGY FOR MINIMUM TIME PROBLEM 
WITH TORQUE CONSTRAINTS ONLY 

Curve A is potential energy 
Parabola B is deformation constraint. 

Potential Energy 

40 00 _ i 0.001 



0.00 0.10 0.20 0.30 0.40 030 0.60 0.70 0.10 0.90 1.00 

Normalized Time 
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POTENTIAL ENERGY FOR MINIMUM-TIME PROBLEM 
WITH TORQUE AND POTENTIAL ENERGY CONSTRAINTS 


Curve A is potential energy 
Parabola B is deformation constraint. 


Potential Energy 
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OPTIMAL CONTROL FOR MINIMUM-TIME PROBLEM 
WITH TORQUE AND POTENTIAL ENERGY CONSTRAINTS 


Note: The optimal final time is 0.8177 seconds, an increase of only 3.7 

percent over the solution of P^. 


Control 
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